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We investigate the polarized parton distribution functions (PDFs) and their uncertainties by using 
the world data on the spin asymmetry A\. The uncertainties of the polarized PDFs are estimated by 
the Hessian method. The up and down valence-quark distributions are determined well. However, 
the antiquark distributions have large uncertainties at this stage, and it is particularly difficult to 
fix the gluon distribution. The \ 2 analysis produces a positively polarized gluon distribution, but 
even Ag(x) = could be allowed according to our uncertainty estimation. In comparison with the 
previous AAC (Asymmetry Analysis Collaboration) parameterization in 2000, accurate SLAC-E155 
proton data are added to the analysis. We find that the E155 data improve the determination of the 
polarized PDFs, especially the polarized antiquark distributions. In addition, the gluon-distribution 
uncertainties are reduced due to the correlation with the antiquark distributions. We also show the 
global analysis results with the condition Ag(x) — at the initial scale, Q 2 = 1 GeV 2 , for clarifying 
the error correlation effects with the gluon distribution. 

PACS numbers: 13.60.Hb,13.88.+e 



I. INTRODUCTION 

The spin structure of the nucleon has been investigated 
extensively since the discovery of the EMC spin effect 
I]. Despite a naive expectation that the nucleon spin 
is carried by quarks, the experimental result indicated 
that only a small fraction is carried by the quarks. In 
order to determine this quark spin content and internal 
spin structure, the determination of the first moments 
of the polarized parton distribution functions (PDFs) is 
necessary. Furthermore, the x and Q 2 dependence of 
these functions is crucial in the calculation of the po- 
larized cross sections. The functions should be deter- 
mined from global analyses of polarized experimental 
data. Such analyses have been made by several groups 
HSSSSSHSHIHGim. Now, there are avail- 
able data for the spin asymmetry A\ or the structure 
function 171 by polarized deep inelastic scattering (DIS) 
experiments [H lill ITU Hrj \u\ - These data are valuable 
especially for fixing the the polarized valence quark dis- 
tributions Au v {x) and Ad v (x). The polarized antiquark 
distributions are still not well determined from the data. 
In particular, their flavor dependence is not shown re- 
liably at this stage. Furthermore, the polarized gluon 
distribution cannot be fixed from the polarized DIS data 
although there are some constraints. 

The above statements describe the current status of 



*E-mail: mhirai@rarfaxp.riken.jp 
tE-mail: kumanos@cc.saga-u.ac.jp 
t E-mail: saito@nh.scphys.kyoto-u.ac.jp 




global analyses qualitatively well; however, we have been 
longing for more quantitative discussion on the uncer- 
tainties in the PDFs. Recently, uncertainty estimation 
methods have been developed for the unpolarized PDFs. 
Mathematical formulations of PDF uncertainties are pro- 
posed in Ref s llsLlial . Practical methods are also devel- 
oped [20, ml 122. 123 and the y a re actually used in recent 
unpolarized analyses |24ll25ll2q . Thanks to a large num- 
ber of experimental data points with excellent precision 
and wide kinematical coverage, the unpolarized PDFs are 
well determined from small x to large x with a reasonable 
precision so that hadronic cross sections can be calcu- 
lated to a few percent accuracy 0> HE • I n t ne same 
way, the uncertainties of the polarized PDFs have been 
investigated 0, 0. The polarized gluon distribution 
has large uncertainties due to a lack of data which are 
sensitive to the distribution. However, because polarized 
experimental projects are going on, we expect to have 
better determination of Ag in the near future. The de- 
termination of the polarized PDF uncertainties enables 
the uncertainty estimation of physical observables such 
as scattering cross sections and spin asymmetries. The 
uncertainty estimation of the polarized PDFs is valuable 
for understanding the present situation and for clarifying 
the importance of future experiments. 

There are three major purposes in this paper. First, 
the uncertainties of the polarized PDFs are investigated 
for the AAC (Asymmetry Analysis Collaboration) pa- 
rameterization [jj. Although there are uncertainty esti- 
mations for the polarized PDFs [13, IT^ . the results could 
depend on the parameterized functional form and the de- 
tails of the uncertainty estimation method. Therefore, it 



is important to estimate the PDF uncertainties indepen- 
dently. In particular, we discuss the large uncertainties of 
the polarized gluon distribution. Second, we investigate 
the impact of precise SLAC-E155 proton data, which are 
not included in the previous AAC analysis, on the PDF 
uncertainties, especially on those of the polarized anti- 
quark and gluon distributions. Third, error correlation 
is investigated by a global analysis with Ag(x) = at 
the initial Q 2 point. We compare its PDF uncertainties 
with those of the Ag(x) ^ analysis in order to show 
error correlation effects. 

This paper is organized as follows. In Sec. [H] we 
describe the framework of our parameterization for the 
polarized PDFs. The Hessian method is explained in 
Sec. Mil as an uncertainty estimation method for the 
PDFs. In Sec. IIVI \ 2 analysis results are shown with the 
polarized PDF uncertainties. First, they are compared 
with the experimental data for the spin asymmetry A\. 
Second, obtained polarized PDFs are shown in compar- 
ison with the distributions of the previous AAC version 
and other parameterization studies. Third, effects of the 
SLAC-E155 data are explained, and the correlation be- 
tween the antiquark and gluon distribution uncertainties 
is discussed. The results are summarized in Sec. [V] 

II. PARAMETERIZATION OF POLARIZED 
PARTON DISTRIBUTION FUNCTIONS 

The major source of information on the polarized 
PDFs has been polarized electron and muon DIS exper- 
iments. The polarized PDFs are determined by compar- 
ing theoretical functions with the asymmetry A\ data of 
the polarized DIS experiments 0, El, EE EE G3 The 
variable Q 2 is given by Q 2 = —q 2 with the momen- 
tum transfer q, and the scaling variable x is defined by 
x = Q 2 /(2p ■ q) with the nucleon momentum p. We se- 
lected the data with Q 2 > 1 GeV 2 so that perturbative 
QCD could be applied relatively safely. Then, the to- 
tal number of available data is 399, and they cover the 
kinematical region, 0.004 < x < 0.75 and 1 < Q 2 < 72 
GeV 2 . 

The spin asymmetry A\ is expressed in terms of the 
polarized structure function g\, the unpolarized struc- 
ture function F2, and the longitudinal-transverse struc- 
ture function ratio R: 

At(x, Q 2 ) = ^ 2 x [1 + R(x, Q 2 )} . (1) 

The SLAC analysis results |23 are used for the function 
R(x,Q 2 ). The structure function i^fx, Q 2 ) is expressed 
by unpolarized PDFs: 

n f r 

F 2 (x, Q 2 ) = J2 e t x \ C a( x i a s) ® h(%, Q 2 ) + ft (a, Q 2 )} 
+ C g {x,a s )®g{x,Q 2 )\. (2) 



Here, q(x,Q 2 ), q(x,Q 2 ), and g(x,Q 2 ) are the quark, 
antiquark, and gluon distributions, respectively, and 
C q (x,a s ) and C g (x,a s ) are coefficient functions. The 
symbol (g) denotes the convolution integral: 

f(x)® 9 (x)= f *Lf(*)g(y). (3) 

Jx V \VJ 

In the same way, the structure function g\(x,Q 2 ) is ex- 
pressed as 

9i(x,Q 2 ) = \ll e2 \ AC q {x,a s )® [Aq t (x,Q 2 ) 
i=l 

+ Aq t {x,Q 2 )] + AC g {x,a s ) ® A<7(x,g 2 )|, (4) 

where Aq(x,Q 2 ), Aq(x,Q 2 ), and Ag(x,Q 2 ) are the po- 
larized quark, antiquark, and gluon distributions, respec- 
tively. The function Aq is defined by Aq = q^ — q- 1 . which 
indicates the difference between the distribution of quark 
with hclicity parallel to that of parent nucleon and the 
one with hclicity anti-parallel. The functions AC q (x,a s ) 
and AC g (x, a s ) are the polarized coefficient functions. 

As the polarized PDF at the initial scale Qg, we choose 
the functional form: 

Af(x) = [8x v - k{x v - x»)]f{x) , (5) 

where f[x) is the corresponding unpolarized PDF. This 
form is taken for imposing the positivity condition and for 
reducing correlations among the parameters. Optimized 
PDFs are four distributions: Au v (x), Ad v (x), Aq(x), 
and Ag(x), which are defined at Qq by Eq. J5J. The 8, 
v, k, and fj, are free parameters, which are determined by 
a x 2 analysis. 

In principle, the separation of these quark distributions 
can be arbitrarily chosen. For example, alternative choice 
would be Au + (x), Ad + (x) and As + (x), where 

Au + (x) = Au(x) + Au(x) : (6) 

and similar expressions for Ad + (x) and As + (x). Here, 
these quark distributions can be related to ours as 

Au + {x) = Au v (x) +2Aq(x) , 

Ad+{x) = Ad v (x) + 2Aq(x), (7) 

As + {x) = 2Aq(x). 

Here Aq(x) can simply be interpreted as an average of 
polarized sea-quark distributions. 

Practical difference appears, however, when we apply 
constraints on the first moments of the quark distribu- 
tions from the axial coupling constants of octet baryons. 
By denoting the first moments by L Af{x)dx — AF, 
these moments should be connected to 

AU+-AD+ = 1.267 ±0.011, 
AU+ + AD+ - 2AS+ = 0.585 ±0.025. (8) 
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These relations can be rewritten by using our separation 
of quark distributions, 



where 



AU V - AD V = 1.267 ± 0.011 - 2A 2 , 
AU V + AD V = 0.585 ± 0.025 - 4A 3 . 



A 2 ee AU — AD , 

A 3 . ATJ + A °-AS. 



(9) 



(10) 



Since there is no experimental guidance on the size of the 
first moments of the flavor asymmetric distributions (A 2 
and A3), we continue to neglect them in this paper as was 
done in our previous one, although we are also preparing 
a new calculation with those breaking parameters. This 
point will be discussed later in this section. 

The polarized distributions are numerically evolved 
to the Q 2 points of experimental data by the DGLAP 
(Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) evolution 
equations ,2Q] in order to calculate \ 2 - The total x 2 is 
defined by 



[Af^(x,Q 2 )-A<j^(x,Q 2 )} 2 
[AAff (x, Q 2 )f 



E 



(11) 



where A^4f ata is the experimental error including 
both systematic and statistical errors: (AA^ ata ) 2 = 
(AAf at ) 2 + (A^f st ) 2 . The total x 2 is minimized by 
the CERN subroutine MINUIT |30j . 

Here, the systematic errors A^ yst are fully included. 
It would be ideal to include uncorrelated and correlated 
systematic errors separately so that we can perform fully 
consistent uncertainty analysis. The issue of the corre- 
lated errors in the global analysis is indeed inve stig ated 
in the recent unpolarized PDF parametrizations |24ll25|. 
In the polarized PDF analysis, however, these errors are 
not listed separately in papers and it is very difficult to 
access such information. Because of this incompleteness, 
our analyses overestimate the uncertainty in the PDFs. 

In order to obtain a rough picture of the effects of sys- 
tematic uncertainties, we also performed x 2 calculation 
with the statistical uncertainties only, which resulted in 
^20% increase in the \ 2 - I n average, this increase in x 2 
corresponds to the fact that quadratic sum of statistical 
error and full systematic error is larger by ~10%. By 
looking at the individual data points, the increase in the 
errors ranges from a few percent to 50%, which is not in- 
consistent with the average picture. From this exercise, 
we can conclude that our uncertainties of the PDFs are 
overestimated, but only by ~10% in average, although 
we now emphasize the needs of separated systematic er- 
rors so that we can perform fully consistent uncertainty 
analyses in future. 

In comparison with the previous AAC analysis jg , the 
SLAC-E155 proton data are added to the \ 2 analysis. 
In order to demonstrate the impact of this new data set, 



we used the same configuration with the previous analy- 
sis. The renormalization scheme is the MS scheme in the 
next-to-leading order (NLO). The initial scale is set at 
Ql = 1 GeV 2 . The number of flavor is fixed at N f = 3 
and heavy flavor contributions are neglected. 

We use the GRV98 NLO parameter set as the unpo- 
larized PDFs [3l|. Even if other unpolarized PDFs, for 
example CTEQ6 |24|, are used, the results do not change 
significantly. In particular, the polarized gluon distribu- 
tion is modified; however, it is well within the uncertain- 
ties for Ag(x) in Sec. If V Bl We should also mention that 
unpolarized PDF uncertainties are not included in our 
analysis for estimating the polarized PDF uncertainties. 
The value of Aqcd is taken from the GRV unpolarized 
PDF analysis: Aq^ 3 = 299 MeV. Unpolarized experi- 
mental data are generally more accurate than polarized 
ones, so that Aqcd had better be determined by the 
unpolarized analysis. 

In the following, we discuss two important constraints, 
the positivity and flavor-symmetric conditions, on the 
polarized PDF determination. 



A. The positivity condition 

The positivity condition means that the magnitude of 
a polarized cross section should be smaller than the cor- 
responding unpolarized one: \Aa\ < a. In the leading or- 
der (LO), this relation indicates |A/(x)| < f(x) because 
probabilistic interpretation can be applied for the parton 
distributions. However, the condition is not strictly satis- 
fied in the NLO because of higher-order corrections [12 • 
The correction due to the coefficient functions is small 
in the limit x —* 1. Therefore, the positivity condition 
|A/(x)| < f{x) could be used practically for constraining 
the polarized PDFs at large x. 

If this condition is not taken into account in 
the x 2 analysis, it tends to be violated significantly 
in the polarized antiquark and gluon distributions: 
[\Af(x)\/ f(x)] x -+i 3> 1. It could lead to an unphysi- 
cal cross section: |Ac| > a. This unfavorable behavior 
comes from the lack of accurate experimental data in the 
large- a; region. Furthermore, experimental data indicate 
that the spin asymmetry A\ increases monotonically as a 
function of x in the large- a; region. It easily leads to un- 
physical results without the positivity constraint. There- 
fore, the positivity condition for the polarized PDFs is 
practically a useful constraint for avoiding the unphysi- 
cal results, and we decided to impose this condition in 
the x 2 analysis. The condition restricts the range of the 
parameter 6 which controls the large- a; behavior of the 
polarized PDFs: -1 < S < 1. 



B. Flavor symmetric antiquark distributions 

It is now known that unpolarized antiquark distribu- 
tions, u, d, and s, are different from lepton scattering 
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and Drell-Yan experiments [3j|. There are model ex- 
planations, for example, by meson clouds, chiral soliton, 
and Pauli exclusion. These models are extended to the 
polarized antiquark distributions. There are available ex- 
perimental data which may indicate the polarized flavor 
dependence [1,13113; however, they are not accurate 
enough to find the difference between Au and Ad. There- 
fore, there is no reliable data for fixing the difference 
between the polarized antiquark distributions, and the 
determination of the flavor asymmetric distributions still 
depend on separation models @, 0, ^| . 

Even in the alternative quark separation, Au + (x), 
Ad + (x) and As + (x), we cannot address these specific 
questions. In future, the flavor dependence of the po- 
larized antiquark distributions will be investigated, for 
example, at RHIC by W production . If these data 
become available, it makes sense to introduce the flavor 
dependent parameters into the x 2 analysis. 

With the flavor symmetric assumption, i.e. A 2 = 
A3 = 0, the first moments of Au v and Ad v are fixed 
from Eq. ©: AU V = 0.926 and AD V = -0.341. 



III. PDF UNCERTAINTY 

The uncertainties of the parton distributions are esti- 
mated by the Hessian method, which has been used as a 
general statistical method for estimating errors. The un- 
certainties come from measurement errors in the global 
PDF analysis. 

The parameters, e.g. 8, v 1 n, and \x for each distri- 
bution in our analysis, are denoted (i=l, 2, • • •, N), 
where N is the number of optimized parameters. Ex- 
panding x 2 around the minimum point a and keeping 
the leading quadratic terms, we have 

A X 2 =X 2 (a + 6a) - X 2 («) - J2 H v Sa * Sa 3 > ( 12 ) 

where the Hessian Hij is the second derivative matrix 
in the function x 2 ( a )- The first derivative terms do not 
exist because they vanish at the minimum point. 

For estimating the PDF uncertainty, we should eval- 
uate a confidence region of the multivariate normal dis- 
tribution for optimized-parameter errors. Equation l|12|) 
indicates the local behavior around a, and the confidence 
region could be identified by an ellipsoid which is defined 
by Ax 2 - Assuming parabolic curves for the function 
X 2 (a), we can vary the ellipsoid to an arbitrary confi- 
dence level by choosing Ax 2 • In our estimation, the A^ 2 
value is obtained by the following procedure. 

The confidence level P for the normal distribution with 
N degrees of freedom is identified with the one for the 
X 2 distribution: 

P -F*m (^"^ (13) 



where T(m) is the Gamma function. It can be chosen 
to be the probability of one-c-error range of the normal 
distribution (P — 0.6826) for our study in order to com- 
pare with the experimental errors. In the case of one 
parameter (N — 1), we obviously have Ax 2 = 1 from 
Eq. (|13[l . Therefore, Ax 2 — 1 could be simply used 
for calculating the uncertainty if the parameter number 
is one. However, the parton distributions are provided 
with many parameters, so that the Ax 2 value should 
be re-evaluated. For example, the parameter number is 
eleven (N = 11) in our new polarized PDF analysis, and 
it leads to Ax 2 = 12.647. 

The uncertainty of a parton distribution F(x, a) with 
respect to the optimized parameters a is then calculated 
by using the Hessian matrices and assuming linear error 
propagation: 

,(14) 

For the PDF uncertainty estimation, we can analytically 
calculate the gradient terms dF(x,a)/da,i at the initial 
scale Qq. For the estimation at arbitrary Q 2 , each gradi- 
ent term is evolved by the DGLAP evolution kernel, and 
then the PDF uncertainties SAf(x,Q 2 ) are calculated. 
The uncertainties of the polarized structure functions g\ 1 
g™, and gf are calculated by the convolution integrals of 
the PDF gradient terms with the coefficient functions. 

IV. RESULTS AND DISCUSSION 

We report our analysis results. Because one of the ma- 
jor purposes is to show the PDF uncertainties, we ana- 
lyzed the data only for the NLO set. In our new analysis, 
the best fit is obtained with X 2 (Jd.o.f.) = 346.5 (0.893). 



TABLE I: Numbers of the Ai data with Q 2 > 1 GeV 2 and x 2 
values are listed. The notations p, n, and d indicate proton, 
neutron, and deuteron, respectively. 


data set 


No. of data 


x 2 


EMC (p) 


10 


4.5 


SMC (p) 


59 


54.0 


E130 (p) 


8 


4.9 


E143 (p) 


81 


61.1 


E155 (p) 


24 


24.2 


HERMES (p) 


19 


19.1 


SMC (d) 


65 


56.5 


E143 (d) 


81 


93.6 


E155 (d) 


24 


20.3 


E142 (n) 


8 


2.6 


E154 (n) 


11 


3.6 


HERMES (n) 


9 


2.2 


total 


399 


346.5 
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TABLE II: Obtained parameters by the NLO analysis (AAC03). 



distribution S v k fi 

Au v 0.975 ± 0.099 0.000 (fixed) 0.601 1.095 ± 0.266 

Ad v -1.000 ± 0.377 0.000 (fixed) -0.721 1.318 ± 0.466 

Aq 1.000 ± 0.994 1.014 ± 0.182 - 90.96 ± 13.57 1.000 (fixed) 

Ag -1.000 ± 3.959 2.248 ± 1.089 254.2 ± 180.7 2.217 ± 2.172 



The x 2 contributions from all the used data sets are listed 
in Table [J and optimized parameters are summarized in 
Table CU 

In the previous AAC version (AAC00) [g, we found 
that an antiquark parameter /i^, which determines the 
functional behavior of Ag at small x, cannot be fixed 
from the data because of the lack of small- a; data. There- 
fore, we fixed the parameter at [iq = 1 in our new analy- 
sis. Other four parameters are also fixed. The parameter 
v Uv and i>dv tended to stop at the positivity limit, so 
that these parameters are finally fixed. The parameters 
K Uq] and are determined by the first moments Au v 
and Ad v from semileptonic data with the assumption for 
the flavor symmetric antiquark distributions. The dif- 
ference from the AAC00 NLO-2 analysis is the addition 
of the SLAC-E155 data. In order to discuss the influ- 
ence of such accurate experimental data on the polarized 
PDF analysis, new analysis results are compared to the 
AAC00-NLO-2 in following subsections. The total num- 
ber of the optimized parameters is eleven, so that the 
uncertainty is estimated by A\ 2 — 12.647 as explained 
in Sec. HTH 



A. Spin asymmetries 

We discuss the results for the spin asymmetries and 
their uncertainties. In addition to the data used for the 
AAC00 analysis, the SLAC-E155 proton target data are 
included. The E155 data cover the region, 0.015 < x < 
0.75 and 1.22 < Q 2 < 34.72 GeV 2 . Calculated spin 
asymmetries are shown in Fig. ^ and they are com- 
pared with the previous results (AAC00 NLO-2). The 
solid curves and shaded areas show the spin asymme- 
tries A\, A™, and Af and their uncertainties of the new 
results (AAC03), respectively. The dashed and dashed- 
dot curves indicate those of the AAC00 NLO-2. The A x 
uncertainties are calculated by using the estimated un- 
certainties for the obtained polarized PDFs. 

The spin asymmetries are slightly modified especially 
in the region 0.02 < x < 0.3. It is noteworthy that 
the E155 proton data also affect the spin asymmetries 
of the neutron and deuteron since Au v (x) and Aq(x) 
are modified. Although the asymmetry modifications are 
rather small, the uncertainties are significantly modified. 
Comparing the shaded areas with the dashed-dot curves, 
we find that the addition of the E155 data reduces the 
Ai uncertainties. The uncertainties in the region x < 0.6 
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FIG. 1: Calculated spin asymmetries and their uncertainties 
are shown at Q 2 = 5 GeV 2 . The solid curves and shaded 
areas indicate the spin asymmetries and their uncertainties 
of the new results (AAC03), respectively. The dashed and 
dashed-dot curves show those of the previous results (AAC00 
NLO-2). 



are reduced directly by the El 55 data. In addition, the 
data indirectly contribute to the uncertainty reduction in 
the large-x region, where precise data are not available, 
through the x-dependent PDF form. 

Next, the differences between the A\ data and the the- 
oretical asymmetries, namely Af ata — A^ ACQ3 , are shown 
in Fig. |21 The shaded areas indicate the uncertainties of 
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FIG. 2: Comparison of the AAC03 spin asymmetries with experimental data. The ordinates indicate the differences between 
experimental data and theoretical values (Af ata — A^ AC03 ). The error bars indicate the errors obtained by the quadratic 
summations of the statistical and systematic errors. The shaded areas show the uncertainties at Q 2 = 5 GeV 2 . 



the AAC03 analysis. The error bars indicate the ones 
obtained by the quadratic summations of the statisti- 
cal and systematic errors. We find that the uncertain- 
ties are roughly equal to the errors of the experimental 
data. The spin asymmetries are constrained in the re- 
gion x < 0.1; however, they still have rather large un- 
certainties in the larger- a; region. It is obvious that the 
large-x asymmetries are not determined well from the 
present data. These uncertainties at large x imply that 
the positivity of the spin asymmetry is not necessarily 
guaranteed in the x 2 analysis unless the positivity con- 
dition is enforced. The shaded areas spread out in the 
region 0.1 < x < 0.6 due to large errors of the E143 and 
SMC data for the proton. The numbers of these data are 
larger than those of HERMES and E155 experiments, so 
that their overall \ 2 contributions are larger and the ac- 
curate E155 data cannot contribute much in this region. 
The situation is similar for the deuteron uncertainties. 
The neutron uncertainties are still large because the used 
3 He target data are not accurate enough in comparison 
with the proton and deuteron data. We expect that the 
neutron uncertainties could be improved by precise JLab 
measurements l37l. 



B. Polarized parton distribution functions 

We discuss obtained polarized PDFs and their uncer- 
tainties. In Fig. El the polarized PDFs of the AAC03 
analysis are compared with those of the with the AAC00 
NLO-2 at Q 2 = 1 GeV 2 . The distributions Au v , Aq, and 
Ag are slightly modified, but Ad v is scarcely changed. 
This is because Au v is the dominant component of the 
proton structure function g\ and because the whole sea- 
quark contribution is in general larger than the Ad v con- 
tribution. Therefore, the larger components Au v and Aq 
are mainly affected by the added precise E155 data. 

Next, PDF uncertainties are shown at Q 2 = 1 GeV 2 
in Fig. 0] The shaded areas are the uncertainties of the 



AAC03 analysis, and the dashed curves indicate those 
for the AAC00. The valence-quark distributions are suf- 
ficiently constrained by the polarized DIS data. How- 
ever, we find rather large uncertainties in the region 
0.1 < x < 0.6, which corresponds to the region of the 
large A\ uncertainties in Fig. |21 It indicates the necessity 
of accurate data in this kinematical region for better de- 
termination of the polarized valence-quark distributions. 
In particular, accurate 3 He data are useful for reducing 
the Ad v uncertainties because the Ad v contribution to g™ 
is almost the same as the Au v (g™ oc AAd v + Au v + 12Aq) , 
whereas the contribution is small in g\. On the other 
hand, the uncertainties of the antiquark and gluon dis- 
tributions are still large. The huge gluon uncertainties 
indicate that the present data cannot rule out the pos- 
sibility of Ag(x) — and negative gluon polarization, 
although the obtained gluon distribution is positive. 




FIG. 3: Obtained polarized parton distributions at Q 2 — 
1 GeV 2 . The solid curves indicate the new AAC03 results, 
and the dashed curves are taken from the previous analysis 
(AAC00 NLO-2). 
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FIG. 4: Polarized PDF uncertainties are shown at Q 2 = 1 
GeV 2 . The solid curves and shaded areas are the polarized 
PDFs and their uncertainties of the new AAC03 results, and 
the dashed curves are the uncertainties of the AACOO results. 



As shown in Fig. all the PDF uncertainties are 
significantly reduced in the AAC03 analysis in compari- 
son with the AACOO because the accurate E155 data are 
added to the data set. In particular, the Ad v uncertain- 
ties are reduced although the Ad v distribution stays al- 
most the same. In addition, the uncertainties of the anti- 
quark and gluon distributions are significantly improved. 
The antiquark uncertainty reduction is directly due to 
the E155 data. However, it is difficult to understand 
that the significant reduction of the gluon uncertainties 
is due to the added new data. This is because that the 
gluon distribution indirectly contributes as a higher order 
correction with the coefficient function, and this contri- 
bution is less than quark contributions. The huge gluon 
uncertainties explicitly indicate the difficulty of fixing the 
gluon distribution from the DIS experimental data. 

We find that the gluon uncertainty reduction is caused 
by an error correlation. The non-diagonal part of the 
Hessian indicates a strong correlation between the polar- 
ized antiquark and gluon distributions. The correlation 
affects the determination of these distributions. We dis- 
cuss the details of the uncertainty improvement due to 
the correlation in Sec. IIV El 



C. Quark spin content 

We show the first moments of the AAC03 parameter- 
ization at Q 2 = 1 GeV 2 in Table ITTT1 and they are com- 
pared with those of the AACOO NLO-2 set. The first 
moments of the up- and down-valence quark distribu- 
tions are fixed in both analyses. The first moments indi- 
cate that quarks carry about 20% of the parent nucleon 
spin, and gluons carry a large positive fraction of the nu- 
cleon spin. Their uncertainties are significantly reduced 
by the added E155 data; however, the present data are 
not enough to obtain accurate values, especially for the 



TABLE III: The first moments of the obtained polarized 
PDFs at Q 2 = 1 GeV 2 . The AAC03 analysis results are 
compared to those of the previous results (AACOO NLO-2). 
The AE is the quark spin content. 





Aq 


Ag 


AE 


AAC03 


-0.062 ± 0.023 


0.499 ± 1.266 


0.213 ± 0.138 


AACOO 


-0.057 ± 0.037 


0.533 ± 1.931 


0.241 ± 0.225 



gluon first moment. 

The uncertainty of the spin content AS strongly de- 
pends on the antiquark uncertainty because it is given by 
AEjv /= 3 = Au v + Ad v + 6Aq. The first moments of the 
valence-quark distributions are fixed, so that the AS un- 
certainty is equal to six times the Aq uncertainty, which 
could be large due to the uncertainty of the distribution 
Aq(x) in the small- x region. It suggests that the extrap- 
olation into the smaller- a; region should be ambiguous 
in calculating the integral over x. We expect that accu- 
rate polarized antiquark distributions will be measured 
in future, then the quark spin content issue will become 
clear. 



D. Comparison with other parameterizations 

The AAC03 analysis results are compared with other 
parameterizations at Q 2 = 1 GeV 2 in Fig. [SJ For com- 
parison, we choose three sets of the polarized PDFs in 
the NLO: GRSV01 (standard scenario) 0, BB (ISET=3) 
[P3| , and LSS (MS scheme) ^3 ■ These parameterizations 
used basically the same experimental data set of the po- 
larized DIS, but they choose averaged data tables over x 
and Q 2 , whereas full tables are used in the A AC analysis. 



Q 2 = 1 GeV 2 




AAC03 


■ ■ ■ GRSV 




- - BB 




-■■ LSS 
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i- 




3 _ 
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4: — 
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FIG. 5: The AAC03 PDFs at Q 2 = 1 GeV 2 are compared 
with the ones for other parameterizations by GRSV01 (stan- 
dard scenario) 0, BB (ISET=3) O, and LSS (MS scheme) 
|T^| . The shaded areas are the uncertainties of the AAC03 
analysis. 



The flavor symmetric antiquark distributions are used in 
all these parameterizations. 

Because the first moments are fixed in the same way, 
the variations are small in the polarized valence-quark 
distributions among the parameterizations in Fig. [S] 
There are slight variations in the antiquark distributions, 
and the gluon distributions differ significantly among the 
analysis groups. However, we find that all the parame- 
terizations are consistent each other because the distri- 
butions are mostly within the estimated error bands. 

The BB and LSS groups also investigated the polar- 
ized PDF uncertainties by the Hessian method. How- 
ever, these uncertainties may not be directly compared 
with our uncertainties because the used A\ 2 value is not 
specified in their papers. In addition, uncertainty esti- 
mation methods are slightly different; for example, only 
the statistical errors are used and a relative normaliza- 
tion shift is introduced in the BB analysis. In general, 
the error estimations involve complicated systematic er- 
rors, e.g. functional form, data selection, and higher- 
twist effects, in the global analyses, and they may not be 
estimated numerically. It is difficult to clarify these is- 
sues only by the current polarized DIS experiment data. 
Therefore, we need to investigate such hidden uncertain- 
ties when we improve the quality of the polarized PDFs 
by incorporating future experimental data. 

E. Ag(x) = analysis and error correlation 

In order to understand the reduction of the gluon un- 
certainty in Fig. 0J we discuss the error correlation be- 
tween the polarized antiquark and gluon distributions. 
As the most simple assumption which could be allowed 
within the gluon uncertainties in Fig. 0] we choose 




0.001 0.01 0.1 1 0.001 0.01 0.1 1 

X X 

FIG. 6: The PDF uncertainties of the Ag(x) = analysis 
are compared with those of the Ag(x) ^ analysis (AAC03). 
The solid curves indicate the polarized PDFs of the Ag(x) — 
analysis at Q 2 = 1 GeV 2 , and the shade areas are their 
uncertainties. The dashed curves indicate the uncertainties 
of the Ag(x) 7^ analysis. 



Ag(x) — at the initial scale (Q 2 = 1 GeV 2 ). However, 
one should note that a finite distribution Ag{x) ^ ap- 
pears at larger Q 2 from the singlet Q 2 evolution. Since 
the gluon-distribution parameters are fixed, we can carry 
out an uncorrelated analysis with the gluon distribution. 

In the Ag(x) = analysis, we obtain x 2 (/d.o.f.) — 
355.0 (0.915), which is 2.5% larger than the value for the 
Ag(x) 7^ analysis. Because it is a slight change in the 
X 2 value, it is reasonable that Ag(x) — is allowed in Fig. 
0] if the uncertainties are taken into account. Obtained 
polarized PDFs are shown in Fig. [5] for the Ag(x) = 
analysis. The total number of optimized parameters is 
seven for this analysis, so that the uncertainties are esti- 
mated by A% 2 = 8.180. The calculated uncertainties are 
shown by the shaded areas, and they are compared with 
those of the Ag(x) ^ analysis (AAC03) shown by the 
dashed curves. We find that the antiquark uncertainties 
are significantly reduced. On the other hand, the valence 
quark uncertainties are scarcely changed, which indicates 
that the correlation with the gluon distribution is weak 
even in the Ag(x) =/= analysis. 

It suggests that the antiquark distributions should be 
determined well by the present data if their errors are 
uncorrelated with those of the gluon. However, because 
of the existence of the strong error correlation, the ac- 
tual antiquark uncertainties are increased by the huge 
gluon uncertainties. In future, if the polarized gluon dis- 
tribution is measured accurately, the antiquark uncer- 
tainties also become small due to the strong correlation. 
To the contrary, the gluon uncertainties could be reduced 
by accurate measurements of the antiquark distributions. 
In this way, we find the significant reduction of the an- 
tiquark uncertainties in Fig. is caused by the error 
correlation effects with the polarized gluon distribution. 
Furthermore, it indicates that the gluon uncertainty re- 
duction in Fig. 0] is also due to the correlation effect 
between these distributions. 

From these studies, it becomes clear that accurate de- 
termination of the gluon distribution is important also for 
the determination of the antiquark distributions. At this 
stage, the polarized gluon distribution is not accurately 
determined, and it also makes it difficult to fix the anti- 
quark distributions from the DIS experimental data. In 
this sense, it is important to measure the polarized gluon 
distribution, for example, by direct photon production 
and jet production at RHIC. 



V. SUMMARY 

We have investigated the optimum polarized parton 
distributions by analyzing the polarized DIS data. We 
focused our studies particularly on three aspects, the un- 
certainty estimation of the obtained PDFs, the role of 
the accurate El 55 proton data, and the error correlation 
between the polarized gluon and quark distributions. 

First, the obtained PDF uncertainties indicated that 
the polarized valence-quark distributions are determined 
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well, that the uncertainties of the polarized antiquark 
distributions arc slightly large, and that the gluon uncer- 
tainties are huge. It is obvious that the polarized gluon 
distribution cannot be determined from the present DIS 
data. 

Second, we discussed the role of accurate E155 data in 
the global analysis. Comparing the AACOO and AAC03 
results, we clarified that the E155 data contributed to 
reducing the PDF uncertainties significantly. 

Third, the error correlation is investigated by repeat- 
ing the parametrization analysis with the initial condi- 
tion Ag(x) — 0. In the Ag(x) = analysis, there is no 
error correlation between the gluon and quark distribu- 
tions, and it leads to small uncertainties of the obtained 
antiquark distributions. This fact suggests that precise 
gluon measurements should be valuable for a better de- 
termination of the polarized antiquark distributions. The 
opposite way is also right: Precise polarized quark mea- 



surements should provide constraints for the gluon dis- 
tribution. 

Finally, we mention that the AAC03 polarized PDF 
library is available at the web site ^ij^. The polarized 
PDFs can be calculated numerically at given x and Q 2 
values. 
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